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Abstract 

The potential for improved performance of wave rotor 
cycles through the use of passage height variation is 
examined. A Quasi-one-dimensional CFD code with 
experimentally validated loss models is used to 
determine the flowfield in the wave rotor passages. 
Results indicate that a carefully chosen passage height 
profile can produce substantial performance gains. 
Numerical performance data are presented for a specific 
profile, in a four-port, through-flow cycle design which 
yielded a computed 4.6% increase in design point 
pressure ratio over a comparably sized rotor with 
constant passage height. In a small gas turbine topping 
cycle application, this increased pressure ratio would 
reduce specific fuel consumption to 22% below the un- 
topped engine; a significant improvement over the 
already impressive 18% reductions predicted for the 
constant passage height rotor. The simulation code is 
briefly described. The method used to obtain rotor 
passage height profiles with enhanced performance is 
presented. Design and off-design results are shown 
using two different computational techniques. The 
paper concludes with some recommendations for further 
work. 

Introduction 

Pressure-gain wave rotors represent a promising 
technology for use as high pressure, high temperature 
topping cycles in gas turbine engines. 1 3 These devices 
utilize unsteady gasdynamic waves, rather than the 
airfoils of their conventional turbomachinery 
counterparts, to transfer energy to and from the working 
fluid. Besides significantly improving the performance 
of gas turbine engines which they top, the potential 
advantages of wave rotors include rotor metal 
temperatures substantially below the combustor 
discharge temperature, rotational speeds which are 
approximately one third those of conventional 
turbomachinery for the same application, a wide 
operating range, 4 and relatively simple rotor geometry. 
Recent research efforts have focused largely on four- 
port, through-flow cycles with axially aligned passages 



Figure 1 Example four-port wave rotor. 

of uniform cross section. This design, shown 
schematically in Fig. 1, is attractive in that, aside from 
the difficulty of partial to full annular transition ducts, it 
can be readily integrated into existing gas turbine 
engines as another spool. 5 

The term ‘through-flow’ refers to the general tendency 
for all of the flow entering the wave rotor to completely 
traverse the passage before exiting. In other words, the 
flow comes in one end and out the other. In contrast, 
there are reverse-flow cycles which draw in gas through 
one port and discharge it through another at the same 
end of the rotor. In these cycles, with reference to 
Fig.l, air would enter the rotor through port 1 and be 
sent to the combustor through port 3. Thus, in a 
reverse-flow cycle each gas stream entering the 
passages tends to remain close to one end of the rotor. 1 

Although the predicted performance benefits of wave 
rotors are impressive, there is always room for 
improvement. With reference to Fig. 1, the higher the 
overall stagnation pressure ratio of the device, P 04 /P 01 , 
the greater the performance benefit to the gas turbine 
that it tops. 1 These performance benefits include 
reduced specific fuel consumption (SFC) and increased 
specific power. Indeed, for small turboshaft engines, 
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with today’s turbine inlet temperatures, there is nearly a 
one-to-one correspondence between SFC reduction and 
wave rotor pressure ratio. That is, a 1% increase in 
P 04 /P 01 yields a 1% decrease in SFC. Thus, for many 
envisioned wave rotor topped gas turbine engines there 
is strong motivation for even small wave rotor 
performance enhancement. 

Significant analytical and computational efforts have 
been directed toward optimizing four-port, pressure- 
gain wave rotors in order to maximize their 
performance; however, until recently, these efforts have 
been constrained to passages which have uniform cross 
section. 3, 6 There seem to be several historical reasons 
for this restriction. The first is that there is a greater 
knowledge base of the one-dimensional flowfields that 
constant area cross sections imply. Secondly, there is 
little to suggest that cross sectional area variation 
should be beneficial . 7 Third, and this is clearly related 
to the second reason, there is nothing to suggest a 
particular manner in which the passage height ought to 
be varied in order to achieve greater performance. 
Even if one could incorporate area variation into the 
analyses, what area profile should be chosen? During 
earlier research periods, when computers were not yet 
available and wave rotor cycles were laboriously 
calculated by hand, the answer to this question was 
simply out of the researchers grasp. Calculation of a 
single cycle was literally weeks worth of work and, 
although incorporating area variation into analysis 
techniques such as the Method of Characteristics was 
reasonably straightforward, any additional degree of 
freedom meant just that many more cycles to calculate. 
The task would have been daunting. Today however, 
the calculation of a design point wave rotor cycle takes 
seconds (or less) on a modem computer and it does not 
seem so overwhelming to attempt at least an exploratory 
foray into the potential benefits of area variation. 
Although there is little information suggesting a 
particular optimal profile, it is relatively easy to try 
many variants and thus make some headway. 


Example passage 
height profile 

-<£ 

Constant passage height 

with equivalent design 



Figure 2 Example of passage height variation. 

benefit to be gained. 

Under the assumption of quasi-one-dimensional 
(heretofore referred to as Q1D) flow of an ideal gas, the 
governing equations will be presented. The numerical 
integration method will then be briefly presented; 
however, this presentation will rely heavily on past 
publications of an experimentally validated numerical 
wave rotor simulation code . 4,8,9 The approach used to 
find an area profile which yielded improved 
performance will then be described. Wave diagrams as 
well as on and off-design performance data will then be 
presented and compared with that from a constant cross 
section passage for a four-port, through-flow cycle 
designed to top a small helicopter engine with a mass 
flow rate of approximately 4.8 lbm/s, and an upstream 
compressor pressure ratio of 7.8. The paper will 
conclude with a discussion of the results, and 
recommendations for further work. 

Governing Equations 

The variation in cross sectional area is assumed to be 
gradual for this study. The passage cross section is 
assumed trapezoidal with constant mean width. Thus, 
as illustrated in Fig. 2, the area variation comes from 
changes in the passage height. Under these assumptions 
and that of a perfect gas, the governing equations for a 
wave rotor passage may be written in non-dimensional 
form as follows 


This paper will present the preliminary results from just 
such an activity and, as will be seen, the results are 
quite promising. It will be shown, to the degree that the 
numerical methods are accurate, that pressure ratios up 
to 4 . 6 % greater than those of the four-port pressure 
exchanger with constant cross section are possible. It is 
emphasized that this is indeed a preliminary result 
because this investigation represents only a ‘first cut’ at 
the potential benefits of area variation. The possible 
area profiles were fairly restricted and, although 
substantial gains were obtained, no general 
methodology or even guideline was found. In other 
words, there is still quite possibly a good deal more 


3w 3F(w) =§ 
dx dx' 
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( 3 ) 


The distance, x' has been normalized by the passage 
length, L. The time, x has been normalized by the 
characteristic wave transit time, L/a , where a is the 
reference speed of sound. The pressure, p' and density, 
p' have been normalized by their respective reference 
values and the axial velocity, u' has been normalized by 
a*. The passage height, h' has been normalized by the 
reference value h . The ratio of specific heats is 
denoted by y. 

In the constant cross sectional area formulation, the 
source vector, S(w) accounts for viscous, heat transfer, 
and leakage effects. Details on the exact form of the 
components may be found in Refs. 8 and 9. In the Q1D 
formulation, these effects are modeled in same manner; 
however, there is an additional component to the 
momentum source term giving it the following form: 

S 2 =Sf on +p'- (4) 

dx' 

Numerical Integration 

The scheme used to numerically integrate Eq. (1) with 
constant cross section has been detailed in Refs. 4 and 
9. It is a Lax-Wendroff based scheme which utilizes 
Roe’s approximate Riemann solver to obtain flux 
estimates at the numerical cell faces. 10 ’ 11 Roe’s method 
utilizes the primitive variables p\ p’, and u’ to obtain an 
approximate solution to a one-dimensional Riemann 
problem with an arbitrary left and right state. Area 
variation (quasi-one-dimensional) can easily be 
incorporated into this method by replacing p and p with 
the products pb' and p'h', respectively. The additional 
term in the source vector, Eq. (4) is dealt with in the 
same manner as other source terms, with the spatial 
height derivative being estimated using central 
differencing. 

Test Cases 

The numerical integration scheme was tested using 
three flow scenarios. The first is an inviscid, steady- 
state, converging-diverging nozzle containing a normal 
shock in the diverging section. While this is not a 
scenario that will be found in wave rotor applications, it 
is convenient because there is a theoretically exact 



Figure 3 Computed and exact mach number 
distributions for a steady-state converging-diverging 
nozzle. 

solution with which to compare results. Figure 3 
compares Mach numbers for the computed and exact 
flowfield solution with a passage height profile 
described by h'=1.25+0.25cos(27tx'). The value of y is 
1.4. Numerical solutions using both 50 and 100 
numerical cells are shown and in both cases agreement 
with the exact solution is excellent. 

The second test flow scenario is more relevant because 
it is unsteady. Unfortunately, there is no exact solution 
to this problem. It was possible however, to compare 
the solution from the present computational scheme 
with two other types of numerical solutions. The first is 
the Method of Characteristics. 7,12 The second is a fully 
two-dimensional, unsteady CFD code. 13 The scenario is 
that of a passage with a cross section described by 
h'=1.0+0.8x'. The inviscid fluid in the passage is 
initially at rest with a non-dimensional pressure and 
density equal to unity. At x=0 the right end of the 
passage is suddenly opened to a reservoir at pressure 
p'=0.471. This initiates a strong expansion wave which 
travels leftward, reflects from the wall, and returns to 
the right end. The velocity at the right end of the 
passage is eventually brought to zero by the reflected 
expansion and at this point the right end of the passage 
is closed. Figure 4 shows a comparison between the 
results from the present Q1D approach, the Method of 
Characteristics scheme, and the two-dimensional code. 
Non-dimensional time traces of pressure along the left 
wall and velocity along the right end of the passage, as 
computed by the three methods, are shown. The two- 
dimensional results were taken from the mid-height 
numerical cells. For reference, contours of pressure in 
the x-t space, computed using the present Q1D scheme 
are shown in the center of the figure. The value of y for 
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Table 1 Variable passage height design constraints. 



p’ u ' 


Figure 4 Traces of pressure at the closed end and 
velocity at the open end of a variable area passage 
with a reflected expansion wave initiated at the right 
end, computed using the Method of Characteristics, 
a 2D CFD code, and the present Q1D CFD scheme. 

this test was 1.4. It is evident that all three 
computational methods agree very well. 

The third test scenario investigated was a simple 
resonance test. A pipe of uniform cross section, closed 
at one end, and filled with an inviscid fluid has a 
fundamental acoustic resonance frequency that is one 
half that of a cone of the same length and with the same 
maximum radius. 14 Furthermore, the overtones for the 
constant cross section pipe are odd multiples of the 
fundamental mode, while they are integer multiples for 
the cone. An effective conical shape can be simulated 
in Eq. (1) by using a parabolic height profile. These 
acoustic tests were conducted with the present 
numerical scheme using a sinusoidal excitation pressure 
perturbation at the open ends. The initial pressure in 
the pipe and cone was unity. The peak-to-peak pressure 
perturbation was 0.002. The results are not shown here; 
however, it was found that the present numerical 
scheme successfully demonstrated resonance for both 
the straight pipe and cone, at the appropriate 
frequencies, for both the fundamental and second 
modes. Higher modes were not examined. It was not 
possible to numerically describe a perfect vertex for the 
cone (i.e. having no area). The area of the numerical 
cell at the tip of the cone was 0.5% of the maximum 
area. The numerical simulations were performed with a 
spatial resolution of Ax'=0.01 and a time step of 
Ax=.005. 

The test cases described above do not imply formal 
validation of the computational method; however, the 
results are at least encouraging. 


Mid-Height Passage Width 

0.34 in. 

Reference Passage Height, h 

0.86 in. 

Rotor Length 

6.00 in. 

Mean Rotor Radius 

3.21 in. 

Rotor Rotational Speed 

16,800 rpm 

Inlet Mass Flow Rate 

4.78 lbm/s 

T 04 /T 01 

2.21 

Inlet Stagnation Pressure 

114.2 psia 

Inlet Stagnation Temperature 

1081.7 R 

Mean Ratio of Specific Heats 

1.353 

Constant Passage Height P 04 /P 01 

1.203 


Height Profile Selection 

As mentioned in the introduction, there do not appear to 
be any clear guidelines in the literature for choosing a 
passage height profile which yields enhanced 
performance. The strategy of this investigation was to 
utilize the most geometrically simple profile possible. 
The motivation for such an approach was simply that 
less complex profiles have less parametric variants to 
consider. 

For all of the profiles chosen, the cycle design method 
described in Ref. 3 was used to obtain circumferential 
port locations, port pressures and temperatures, and 
reference passage heights, subject to the fixed design 
constraints shown in Table 1. Briefly this method 
utilizes the numerical simulation described above to 
perform a multi-level iteration, the convergence of 
which produces a cycle with port placement to properly 
initiate or reflect waves (i.e., only at port openings and 
closings), and port boundary conditions to insure that 
the same amount of mass leaves the wave rotor as enters 
in a single cycle. Since it utilizes the CFD code 
described above, the design method accounts for the 
same losses, namely those due to viscosity and gradual 
passage opening. Heat transfer, leakage, and incidence 
losses are not incorporated in the design method; 
however, they are included in general numerical 
simulations. Bleed flows for cooling of downstream 
turbine components and ducts can also be accounted for 
in the design code; however, bleed was not used in this 
study. 

The fixed design constraints shown in Table 1 were 
taken from an optimized constant area wave rotor 
designed for topping a small helicopter engine. 1,5 As a 
result, the optimal performance enhancement profile to 
be shown is specific to this small engine application. 
The reference height listed in Table 1 is also the 
constant passage height design value. The performance 
figure of merit used for the investigation was the overall 
stagnation pressure ratio, po 4 /poi- The free parameters 
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in the design method were the height profile and the 
expansion ratio, e which is the ratio of assumed constant 
static pressure in port 4 of Fig. 1 to the assumed 
uniform pressure in the passage just before it opens to 
port 4. For each height profile tested the expansion 
ratio yielding maximum performance was found by trial 
and error. It was always within 10% of the optimal, 
constant height value of 0.40. 

Profiles Described by a Single Slope 

The first height profiles examined were described by 

straight converging or diverging lines of various slopes. 

It was found that passages which diverged in the 
direction of flow showed decreased performance 
compared to the constant height passage. Passages with 
converging profiles showed slight improvements with 
overall pressure ratio reaching a maximum of 
approximately 1.8% above the constant height value 
when the slope of the passage height, dh’/dx' was -0.37. 

Profiles Described by Two Slopes and a Vertex 
Height profiles described by two slopes and a vertex 
were examined next. It was found that those profiles 
which converged then diverged in the flow direction 
showed significant performance improvement. In 
particular, a profile with non-dimensional converging 
slope, diverging slope, and axial location of the vertex 
of -1.97, 0.37, and 0.39 respectively raised the overall 
pressure ratio to 3.4% above the constant height value. 

Converging-Diverging-Converging Profiles 
The converging-diverging profiles yielded even higher 
performance when a second, much shorter converging 
section was added after the diverging section. The 
profile of this type yielding the highest performance is 
shown in Fig. 2. The profile has been smoothed using a 
curve-fit procedure which removes the discontinuous 
slopes at the vertices and assures that the slopes are zero 
at the ends of the passage. This was the last profile 
examined; however, it is noted that the examination was 
not exhaustive. Clearly, each new degree of freedom in 
the height profile description greatly increases the scope 
of the parametric analysis and a complete parametric 
study of the converging-diverging-converging profiles 
was beyond the scope of this paper. For reference, Fig. 
2 also shows the constant height profile giving the same 
mass flow rate under the constraints listed in Table 1. 

Design Point Performance 

The design point overall pressure ratio obtained using 
the optimal passage height configuration shown in Fig. 
2 was 1.255. This value was obtained with heat transfer 
effects included and an assumed endwall leakage gap of 
0.005 in. which is the same as that used in the constant 
passage height simulation with which the Table 1 


Table 2 Design point port conditions for a wave 
rotor with and without passage height variation 


Port # 
Variable h 


(Constant h) 

Po To 

Mach # 

Mass Flow 

1 

1.00 1.00 

0.43 

1.00 


(1.00) (1.00) 

(0.51) 

(1.00) 

2 

3.56 2.06 

0.36 

1.50 


(3.26) (2.04) 

(0.38) 

(1.54) 

3 

3.32 2.87 

0.31 

1.50 


(2.97) (2.82) 

(0.40) 

(1.54) 

4 

1.26 2.21 

0.60 

1.00 


(1.20) (2.21) 

(.60) 

(1.00) 



Figure 5 Computed design point wave diagram with 
optimal passage height profile, shown using contours 
of temperature. 


baseline pressure ratio of 1.203 was obtained. The 
computed design point wave diagram is shown in Fig. 5 
using contours of temperature. The areas between the 
contour lines have been gray-shaded such that dark 
grays represent relatively cool regions and light grays 
represent hot regions. The use of temperature contours 
allows waves patterns and a limited number of particle 
paths to be seen simultaneously. Figure 5 is 
unremarkable in that there is nothing obvious which 
distinguishes the wave pattern from that of a four-port 
cycle with constant passage height. For reference, Fig. 
5 was computed using 100 numerical cells with a 
corresponding numerical time step of Ax=.003. 

The design point stagnation conditions and Mach 
numbers for each of the four ports shown in Fig. 1 (and 
Fig. 5) are listed in Table 2. For ports in which there is 
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Figure 6 Schematic of the numerical test bed. 

flow into the wave rotor, the stagnation conditions are 
either specified in the case of port 1, or supplied by 
component sub-models in the case of port 3. 8 For 

outflow ports 2 and 4, the stagnation conditions are 
calculated using a constant area mixing calculation to 
account for losses due to non-uniformities in the port 
flowfield (see for example the velocity distribution in 
Fig. 4). The stagnation conditions and mass flow rates 
in Table 2 have been normalized by those from the 
inlet. For reference, the port states for the constant 
passage height rotor design are also listed in Table 2. 
Here, the differences between constant and variable 
passage height are clearer. The compression achieved 
by the variable height cycle, P 02 /P 01 is 9.2% higher than 
the constant height cycle. Furthermore, the stagnation 
pressure drop across the combustor required to balance 
the cycle drops from 8.9% of the port 2 stagnation 
pressure, for the constant height rotor to 6.7% for the 
variable height rotor. These results are consistent with 
the improved overall pressure ratio of the variable 
height cycle, suggesting a higher compression efficiency 
and a better match between compression and expansion 
processes. 


Off-Design Performance 

The off-design performance will be examined in two 
ways. The first will be a presentation of a performance 
map which encompasses the entire potential operating 
regime of the wave rotor in an engine topping 
environment. The second will be through detailed 
views of wave diagrams at several off-design operating 
points. 

Performance Map 

The performance map was generated using a numerical 
‘test bed’ shown schematically in Fig. 6. The wave 
rotor component was simulated using the scheme 
described earlier. The other components shown in the 
figure were simulated using simple, lumped-capacitance 
models. 4 Discrete spatial and temporal steps of Ax=.02 



(^corrected ) design 

Figure 7 Four-port wave rotor performance map 
with and without area variation. 

and At=.005, respectively were used. The wave rotor 
cycle described in this paper was designed to operate in 
a so-called freewheeling mode. This means that there is 
no drive motor to turn the rotor. The small amount of 
torque required is generated by changes in angular 
momentum as fluid in the inlet ducts is turned from the 
duct angles to follow the walls of the rotor. Windage 
and bearing friction are neglected. Inefficiencies 
associated with the flow turning process are modeled as 
incidence losses. For a given operating point, the rotor 
speed changes until the net torque generated by the inlet 
ports is zero. With this type of design, and the perfect 
gas assumption, the steady-state performance of a wave 
rotor may be displayed on a single map consisting of 
families of curves with fixed corrected heat addition 
rates (fuel flow rates) plotted with corrected mass flow 
on the ordinate and overall pressure ratio on the 
abscissa. Such a map is shown in Fig. 7. The map 
shows performance for both the variable and constant 
height wave rotor designs. The constant height results 
are shown with symbols and lines. The variable height 
results are shown with only lines. For reference, the 
corrected mass flow rate is defined as 


m 


m„ = 


Poi A iV S 


RT, 


01 


( 5 ) 


where m is the mass flow rate, p m and T 0 i are the inlet 
stagnation pressure and temperature respectively, Ai is 
the inlet cross sectional area, R is the gas constant, and 
g c is the Newton constant. The corrected heat addition 
rate is defined as 


Qc = 


Q 

Poi A lV^Olgc 


( 6 ) 
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where Q is the heat addition rate. The corrected heat 
addition and corrected mass flow are shown in Fig. 7 as 
fractions of the design point value in order to facilitate a 
more direct comparison between the variable and 
constant passage height cycles. Also shown in Fig. 7 is 
an approximate operating line, in an engine 
environment. The line represents operation from a near 
idle condition to 100% of power. It can be seen in this 
figure that the variable height design out performs the 
constant height design over the normal operating regime 
in a topping cycle environment. Furthermore, it is clear 
that, like the constant height design, the complete range 
of operation for the variable height design is quite large. 
This is an encouraging result, and parenthetically, 
somewhat counterintuitive given that many 
turbomachinery performance increases come at the 
expense of reduced operating range. 

Off-Design Wave Diagrams 

The off-design wave diagrams were made using a 
steady-state, quasi-three-dimensional CFD code 
(heretofore referred to as Q3D). 15 ’ 16 Details of the Q3D 
formulation may be found in Appendix 1. Briefly, the 
term quasi-three-dimensional refers to a curvilinear 
geometry bound by surfaces of revolution composed of 
the inner and outer surfaces of the wave rotor. The flow 
field is restricted to a surface of revolution midway 
between the two bounding surfaces. Variation normal 
to this mean surface is neglected, and the modeled flow 
becomes two-dimensional in the circumferential, and 
meridional directions. Effects of radius and flowpath 
height variation are thus captured; however no radius 
variation was considered in this investigation. 

Unlike the Q ID code the Q3D code has the capability 
of simulating the interactions between flowfields in the 
wave rotor ducts and those on board the rotor itself. 
These interactions are present in any wave rotor and 
thus, the results presented herein are quite general. It 
has been suggested that these interaction may be quite 
important for accurate flowfield predictions, 
particularly in off design scenarios where inflow to the 
rotor may occur in ducts which are designed for 
outflow, and where pressure gradients across the ducts 
may be large, thus invalidating the constant pressure 
boundary assumption of the one-dimensional, and Q1D 
unsteady codes. 15 This section of the paper will 
therefore serve the dual purposes of illustrating some 
aspects of off-design wave rotor operation, and 
assessing the importance of duct/rotor flowfield 
interactions. Results from the Q3D code will be 
compared to the previously described Q1D code. 
Boundary conditions for ports 1, 3, and 4 in the Q3D 
code were provided by running the Q1D code in the 



Figure 8 Q3D simulation of design point wave 
diagram, with optimal passage height profile, 
neglecting gradual opening time, heat transfer, and 
leakage effects; shown using contours of 

temperature. 

simulation mode described by Fig. 6. The port 2 
boundary conditions were determined by matching the 
mass flow through ports 2 and 3. The Q3D code uses 
sub-models similar to those in the quasi-one- 
dimensional code to capture the effects of viscous and 
incidence losses. No sub-models have been 
implemented to assess leakage, heat transfer, or gradual 
passage opening effects. As a result, the Q1D code 
was, for this portion of the investigation, run with these 
effects disabled. Rotor speed changes associated with 
the free-wheeling mode of operation also have not been 
implemented in the Q3D code. Thus, it is run at a 
specified, fixed speed. For reference, a five zone grid 
was used for all Q3D calculations with 101 (axial) by 
201 (circumferential) cells in the rotor zone, 26 by 49 in 
port 1, 26 by 47 in port 2, 26 by 53 in port 3, and 26 by 
55 in port 4. 

Design Point 

The port timing for the optimal passage height profile 
was designed with the effects of gradual opening time 
accounted for. Since the Q3D code does not have these 
effects modeled, and since finite opening time has a 
substantial effect on the speed of waves in a cycle, it 
stands to reason that the ‘design point’ simulation of 
this cycle would in fact be off-design. In a wave 
diagram, this would be evidenced by waves that do not 
coincide with port openings and closings, and by 
spurious waves in the x-t space. That this is the case is 
evident in Fig. 8. This wave diagram, again made with 
contours of temperature, represents the wave rotor 
operating at approximately the same overall temperature 
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Figure 9 Port static pressure distributions, at the 
design point; on either side of the duct/rotor 
interface for Q3D, in the image cell for Q1D. 

ratio, To/Tqi, and ratio of flow through the combustor 
to flow through the inlet as that of Fig. 5. The contours 
have been shaded such that dark lines correspond to low 
temperature and light lines to high temperature. An 
additional weak shock wave can be seen running from 
the opening of port 1 to just after the closing of port 4 
and reflecting back to the left end of the rotor. It can 
also be seen that the shock initiated at the opening of 
port 3 arrives at the right end of the rotor too early and 
reflects from the wall rather than from the leading edge 
of port 2. The reflected shock running from just below 
the port 2 opening to the port 3 closing is also mis- 
timed, arriving late rather than coinciding with the 
stopping expansion wave which is initiated at the port 3 
closing. 


Despite the improper wave timing, Fig. 8 closely 
resembles Fig. 5 and suggests that the Q ID and Q3D 
codes obtain nearly the same solution to the wave rotor 
flowfield near design operating conditions. This is 
further evidenced in Fig. 9 which shows the computed 
static pressure distributions across each of the ports, at 
the rotor face, for the design operating point. For 
outflow ports 2 and 4 this distribution was specified as 
constant in the Q1D code(as long as outflow was 


Table 3 Performance Prediction Comparison 


Code (case) 

Po7poi 

T 04 /T 01 

rr^/mj 

Q1D (near design) 

1.30 

2.21 

1.51 

Q3D (near design) 

1.28 

2.20 

1.60 

Q1D (idle) 

1.12 

1.60 

1.04 

Q3D (idle) 

1.13 

1.67 

1.10 



Figure 10 Q3D simulation of idle point wave 
diagram, with optimal passage height profile, 
neglecting gradual opening time, heat transfer, and 
leakage effects; shown using contours of 
temperature. 

predicted). For the Q3D code, the static pressure was 
specified at the ends of the ports 2 and 4 and the 
pressure at the rotor face was computed as part of the 
solution. The static pressures at the ends of the ports in 
the Q3D calculations were the same values as those 
assumed at the rotor face for the Q1D results. Figure 9 
shows Q3D results just to the left and right of the rotor 
duct interface, and Q1D results in the image cell which 
may be thought of as being infinitesimally within the 
rotor. For the inflow ports 1 and 3, the Q3D code again 
computed the rotor face static pressure distributions as 
part of the flowfield solution, with the boundary 
conditions being supplied by specified stagnation 
conditions at the entrances to the ports. For the Q1D 
results, the port flows are modeled and the rotor face 
pressure distribution is essentially a boundary 
condition 9 . Again, it is evident from Fig. 9 that the two 
codes predict nearly the same flowfield for this 
operating point. 

Idle Point 

Figure 10 shows the wave diagram representative of an 
idle condition of the wave rotor in an engine topping 
environment. The idle point corresponds approximately 
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Figure 11 Port static pressure distributions, at the 
idle point; on either side of the duct/rotor interface 
for Q3D, in the image cell for Q1D. 

to the bottom of the in-engine operating line of Fig. 7. 
The corrected rotor speed at this point is 86% of the 
design value. The overall temperature ratio, pressure 
ratio, and ratio of combustor to inlet mass flows are 
listed in Table 3. This table shows computed results 
from both the Q1D and Q3D codes for comparison. 
Like Fig. 8, it is clear from Fig. 10 that there are mis- 
timed and spurious waves indicative of off-design 
operation; however, given the operating point, it is 
surprising that the timing is as good as it is. The most 
pronounced feature of the idle wave diagram is the gas 
path. It can be seen that a rather large fraction of the 
fresh air which enters port 1 never reaches the 
combustor, but is instead simply compressed, expanded, 
and sent out through port 4. This leads to the large 
temperature gradients visible in the port 4 flow. The 
impact of such gradients on the performance of 
downstream turbine components is unclear, but it is 
worth noting that they are seen, although less severe, at 
operating points well above idle. 


The port static pressure distributions at the rotor face, 
imposed, or computed by the Q1D, Q3D codes are 
shown in Fig. 1 1 . Again it is clear that there is little 



Figure 12 Streamlines from Q3D simulation of 
severe off design operation, with optimal passage 
height profile, neglecting gradual opening time, heat 
transfer, and leakage effects. 

difference between the two codes, even at this off- 
design operating situation. 

Severe Off-Design Point 

Although off-design situations where port flow reversals 
occur are possible, it was found with both the Q1D and 
Q3D codes that these were only encountered in extreme 
operating conditions, very far from the normal operating 
line of Fig. 5. Nevertheless, it is worthwhile to present 
one such operating point simply to illustrate the 
phenomenon. Figure 12 shows the streamline pattern 
from the Q3D simulation of an operating point in which 
the rotor speed and heat addition rate were held to the 
idle value, but the port 4 back pressure was increased 
until the corrected mass flow rate was 72% of the 
design value. With reference to Fig. 5 it is clear that 
such an operating point is far from normal in-engine 
operation. Regions of flow reversal or nearly stagnated 
flow are visible in Fig. 12 at the tops of ports 1, 3, and 
4. Although not a wave diagram per se, the wave 
pattern can be seen in fig. 12 and can be identified as 
mis-timed. In particular, it is seen that the flow 
reversals in both ports 1 and 3 result from the 
impingement of left running shock waves. It is 
interesting to note that the flow reversal regions were 
correctly predicted by the Q1D code, despite the lack of 
port/rotor interaction simulation capability. This may 
be seen in Fig. 13 which shows the axial velocity 
component in each of the ports at the rotor face, as 
computed by the Q1D and Q3D codes. 
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Figure 13 Port axial velocity distributions, at the 
severe off-design point; on either side of the 
duct/rotor interface for Q3D, in the image cell for 
Q1D. 

Discussion 

The reasons that certain passage height profiles lead to 
enhanced performance are not clear. It is possible, of 
course, that the results are artificial, arising from the 
approximate nature of the quasi-one-dimensional 
assumptions; however, this does not seem likely since 
all of the height profiles examined had small values of 
dh/dx (see Fig. 2). Nevertheless, the results should, at 
least, be validated with a higher fidelity numerical code. 



Figure 14 Plot of initial compression ratio, 7tj C 
versus expansion efficiency, T\ e for several height 
profiles. 

the expansion process, and the initial compression of 
the fresh air. Since incidence, leakage, and heat transfer 
losses are neglected in the design algorithm, 
inefficiencies in the expansion and compression 
processes occur through viscous and gradual passage 
opening losses. Further inefficiencies in the expansion 
process are generated by assumed mixing of the non- 
uniform port 4 velocity profile. Figure 14 shows a plot 
of initial pressure ratio, 7tj C versus the adiabatic 
expansion efficiency, rj e for several of the height 
profiles examined. The initial pressure ratio is defined 
as the ratio of mass averaged pressure in the passage 
after the closing of port 1 to the inlet stagnation 
pressure. The expansion efficiency is defined using the 
mixed stagnation conditions of port 4 and the conditions 
in the passage just before it opens to port 4. Both 7ti C 
and il e are calculated in the relative frame of reference, 
where the design code operates. It was found that all of 
these points could be well described by the following 
function 


In the absence of a complete explanation, it is 
worthwhile to at least make some pertinent observations 
about the wave processes. These observations concern 
the so called ‘bottom half of the four-port wave cycle. 
With reference to Fig. 5, this is the region which begins 
at the bottom of the x-t space and extend to somewhere 
between ports 1 and 3. In this region hot gas is 
expanded and sent out of port 4 while cold, fresh air is 
drawn in through port 1 and initially compressed by the 
coalescing shock wave which originates near the closing 
of port 4 and nearly terminates at the closing of port 1 . 
The overall pressure ratio of the wave rotor, P 04 /P 01 is 
essentially determined by these two processes. The 
design objective is to maximize both the efficiency of 


_ relative 
P04 

„ relative 

P 01 


-1.3127t ic n e -0.325 


(7) 


Equation 7 can be used to obtain constant 
preiative/preiative jj nes ^ two Q f w hich are also shown in 

Fig. 14. Moving to the right of these lines implies 
higher performance. Examining Fig. 14, it is interesting 
to note the relative positions of the various profiles. 
Note, for instance that the diverging profile has 
improved expansion efficiency over the constant height 
profile, but the compression ratio is considerably lower. 
Similarly, the optimal converging profile has lower 
efficiency but higher compression ratio resulting in an 
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overall improvement in pf^ lativc /po e / ativc from the 
baseline constant height profile. The optimal 
converging-diverging profile showed improvement over 
the converging profile largely through increased initial 
compression ratio. The optimal converging-diverging- 
converging profile showed only slight improvement 
over the converging-diverging profile; however, this 
was achieved through a reduction in initial compression 
and a large increase in expansion efficiency. 

It is interesting to note as an aside that the vertex of 
converging diverging profiles, and the first vertex of the 
converging-diverging-converging profiles occurred very 
close to the intersection, in the Fig. 5 x-t space, of the 
port 4 to port 1 coalescing shock and the hot/cold 
interface running leftward from the opening of port 1 . 

Clearly, the above discussion does not satisfactorily 
explain the improved performance observed for certain 
passage height profiles. It is hoped however, that it will 
at least provide a helpful starting point from which 
future studies may benefit. It is certainly clear that 
more work is needed and warranted in this area. 

Concluding Remarks 

Variable passage height has been numerically examined 
in this paper as a means toward improving performance 
in four-port, through-flow wave rotors. Although the 
physical mechanisms are not yet clear, the results 
indicate that substantial improvements can be obtained 
with certain profiles. In particular, an overall pressure 
ratio increase from 1.203 for the constant passage 
height baseline to 1.255 for a height profile described 
primarily by three slopes and two vertices was obtained. 
This substantial improvement suggests that a systematic 
optimization procedure, grounded in an understanding 
of the mechanisms behind the improvements is 
warranted and may lead to even more impressive 
results. 

The off-design simulation results clearly demonstrated 
the robust operational capabilities of 4-port, through- 
flow wave rotors in general, and parenthetically, 
demonstrated the capabilities of and compatibility 
between the Q1D and Q3D computational approaches 
to the flowfield solution. 
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Appendix I Q3D Description 
The Q3D code is similar in physical basis to the model 
presented by Lear and Candler. 15 Finite rotor passage 
widths are neglected, and the limiting condition is then 
steady in the stationary reference frame. With this 
assumption made, flow may be obtained as a steady- 
state boundary value problem. This problem was cast 
using inviscid flow equations and integrated by 
numerical techniques. 

Several differences between the code used here and that 
presented by Lear and Candler exist. First, in the 
present code, the effects of passage height profile 
variation, aerodynamic incidence, and viscosity were 
modeled. Second, the upwind fluxes in the present code 
used Roe flux difference splitting instead of the Steger- 
Warming technique. 11 Third, second order spatial 
accuracy was obtained using van Albada limited 
MUSCL differencing. 18 

Other flux limiters were also considered. The Superbee 
limiter was noted in the literature for its ability to 
resolve contact discontinuities. 18 In testing with the 
Q3D code, resolution of contact surfaces appeared 
better with the Superbee limiter than with the van 
Albada limiter. However, the Superbee limiter yielded 
solutions containing nonphysical oscillations and 
appeared to inhibit steady state convergence. Overall, 
cycle performance predictions were only slightly 
affected by the choice of limiter. 

For cases considered here, within each port, there was 
no profile height, or mean radius variation, therefore, 
the modeled flow was strictly two-dimensional. Thus, 
this flow was modeled with the two-dimensional, 
inviscid, Euler equations. Each port and the rotor was 
numerically integrated on a separate structured grid 
zone. Communication across inter-zone interfaces was 
first order accurate. Boundary conditions were applied 
in a first order manner except for the periodic boundary 
which was second order accurate and consistent with 
the interior fluxes. 

Within the rotor, the theoretical limit of infinite passage 
solidity was applied. Therefore, in the frame of 
reference moving with the rotor, the flow was always 
aligned with the rotor passage. As flow crossed the 
interface between the port and the rotor, in the frame of 


reference of the rotor, it was assumed that no work was 
done. In addition to no work in the rotating frame of 
reference, the assumption was made that some loss due 
to aerodynamic incidence could be applied at the 
interface between the port and the rotor. The 
quantitative form of this loss was applied as a 
discontinuous decrease in relative total pressure for 
flow boarding the rotor across the interface. The loss 
models applied here were applied in a similar manner to 
that described by Paxson. 9 It was also necessary to 
specify the circumferential velocity on the port side of 
the rotor to port interface. For results presented here, 
the circumferential velocity was allowed to convect 
with the axial flow. Thus there was no aerodynamic 
deviation for flow leaving the rotor. Coupled with mass 
continuity, this completed conservation conditions 
across the rotor port interface. 
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